Chirality, charge and spin-density wave instabilities of a two-dimensional electron gas 

in the presence of Rashba spin-orbit coupling 
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We show that a result equivalent to Overhauser's famous Hartree-Fock instability theorem can 
be established for the case of a two-dimensional electron gas in the presence of Rashba spin-obit 
coupling. In this case it is the spatially homogeneous paramagnetic chiral ground state that is 
shown to be differentially unstable with respect to a certain class of distortions of the spin-density- 
wave and charge-density-wave type. The result holds for all densities. Basic properties of these 
inhomogeneous states are analyzed. 



I. INTRODUCTION 



Recent interest in the properties of the quasi-two dimensional electron and hole devices in the presence of structural 
(Rashba-Bychkov)^ or intrinsic (Dresselhaus)^ spin orbit has brought to the fore the problem of the interacting 
chiral electron liquid. It is therefore important to revisit several of the fundamental notions of many-body theory 
for this intriguing system. The purpose of present paper is to begin a theoretical exploration of the relevance and 
special properties of a class of spatially non homogeneous spontaneously broken symmetry states of the electron 
liquid in the presence of Rashba spin-orbit coupling in two dimensions. Specifically we will focus our attention on spin 
density and charge density wave type states henceforth referred to for simplicity as SDW and CDW. SDW and CDW 
states, originally conceived by A. W. Overhauser^^ are generally stabilized by the electron-electron interaction and 
are characterized by spatial oscillations of the spin density, the charge density, or both. In the absence of spin-orbit 
coupling, one can begin to describe SDW and CDW states by simply considering the electron number density for both 
' spin projections: 
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n t = ^+Acos[Q-R+^-\ , (1) 



H = ^+Aco S (^-R-^) . (2) 

> . 

In these expressions the wave vector Q spans the Fermi surface, i.e. does satisfy the condition \Q\ ~ 2kp. A CDW 
corresponds to (f> = 0, while a SDW state obtains for <j> = n. Mixed state are also possible: one such state is beautifully 
realized in chromium^— As we will show, in the presence of linear Rashba spin orbit the corresponding distorted 
states are characterized by a more complex spatial dependence of the number density, the spin density and, where 
appropriate, a chiral density. As a first step towards establishing the fundamental properties of SDW- and CDW-like 
states in the presence of spin-orbit interaction we present here a generalization of the famous Overhauser's Hartree- 
Fock (HF) instability theorem. The latter represents an important exact result in many-body theory for it establishes 
that, within HF, the homogeneous paramagnetic plane wave state does not represents a minimum of the energy for an 
otherwise uniform electron gas for it can be always variationally bettered by a suitably constructed distorted chiral 
SDW or CDW£ 

The paper is structured as follows: In Section [XT] we discuss the relevant aspects of the theory of a two dimensional 
non-interacting electron gas in the presence of Rashba spin-obit coupling. Section Hill briefly discusses useful notions 
of the electron-electron interaction within Hartree-Fock approximation. Section HVl is dedicated to the actual proof of 
the theorem and contains the main results. Finally the last Section contains the conclusions while a number of useful 
mathematical relations are derived in the two Appendices. 

II. TWO DIMENSIONAL ELECTRON GAS IN THE PRESENCE OF RASHBA SPIN-ORBIT 

In the presence of linear Rashba spin-orbit coupling, the one-particle hamiltonian can be written as follows: 

H = Vap- (a x z) , (3) 
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where I is the unit direction along the z-axis, the motion taking place in the x, y plane. 

The non interacting problem can be readily diagonalized to obtain the energy spectrum and the eigenfunctions: 
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where /i = ± labels a state's chirality and 0£ is the angle spanned by the x-axis and the two dimensional wave vector 
k. A schematic of the lower energy sector of the spectrum is plotted in Fig. [1] 




FIG. 1: Non interacting energy spectrum in the presence of linear Rashba spin-orbit coupling. 

By making use of the states of Eqs. §5§ and ([6]) as a basis for a second quantization representation, the familiar 
fully interacting electron gas hamiltonian reads: 



H= E~ b- St +— !— \" yfoklMrf gt p 

S,fi q,klM Ml>^2,M3,M4 

where the phase factor 'Jf 2 ;? „ is defined as follows: 
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Uncorrelated many body wavefunctions for the system at hand can be represented by Slater determinants con- 
structed by occupying any combinations of the chiral states Eqs. (0 and ([6]). 
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III. HARTREE-FOCK THEORY OF A TWO-DIMENSIONAL ELECTRON LIQUID IN THE 
PRESENCE OF RASHBA SPIN-ORBIT COUPLING 



An accurate description of a realistic electronic system requires that electron-electron interaction be taken into 
account. A first step towards developing such a many-body theory is to investigate the results of a mean-field 
approach. The main idea behind the mean field procedure is to find an effective Hamiltonian which is quadratic in 
the electron creation and annihilation operators and can therefore be easily diagonalized. Within the HF theory, 
the ground state is approximated by a single Slater determinant made out of single particle wavef unctions, which 
in turn, are determined by imposing the requirement that the expectation value of the Hamiltonian over the Slater 
determinant be a minimum 10 . Using these wavefunctions as our basis set, a standard Wick decoupling procedure^ 
allows us to determine the effective HF potential. It can be easily proved that the non-interacting chiral states are 
indeed among the solutions of the corresponding HF equations. In this case, the HF potential is diagonal in wave 
vectors and chiral indices: 



Ov=-^#^E^[l + M-os(^-^)] . (9) 

The corresponding HF eigenvalues are given by: 

h 2 k 2 1 

e k^^- a t lk ~^2Y, v k-A l + ^ VC0 <h- ( t ) -)\ ■ ( 10 ) 

An evaluation of the Fermi energy of the two sub-bands leads to an interesting problem. Since one band will, in 
general, acquire more exchange energy than the other, this may result (in a first iteration) in two different Fermi 
levels. In order to equalize them (for elementary stability reasons), electrons from one subband will have to be moved 
to the other. This is the phenomenon of repopulation. 

The spatially homogeneous chiral states are just one of the possible Hartree-Fock solutions. A detailed analysis 
of the possible solutions corresponding to symmetric occupations in momentum space can be done by systematically 
minimizing the total energy as a function of spin orientation and generalized chirality of the system^ 1 -. More general 
solutions correspond to non-symmetric occupations of the single particle chiral states. The problem has been studied 
and the corresponding very interesting phase diagram has been explore d 11 ! 12 . As we will presently discuss there also 
exists an interesting class of spatially non-homogenous solutions to the problem. 



IV. PROOF OF THE INSTABILITY THEOREM 



We will proceed by showing that it is always possible to lower the energy of the homogeneous paramagnetic chiral 
ground state by introducing a suitable real space distortion which is periodic with wave vector Q = 2kpx. The general 
approach follows that of Fedders and Martin^ and is based on an Ansatz which represents a generalization of that 
given by Giuliani and Vignale for the case of the three dimensional electron gas^ 

Let us consider first the putative HF ground state of our many-body system \<&s)- A complete, and, as we shall 
see convenient, description of this determinantal state can be achieved in terms of the corresponding single-particle 
density matrix elements here given by: 

p a p = ($s|&U/»|*sr) , (11) 

where here a and /3 label the one-particle states which are used to build the Slater determinant. 

Now, within the space of Slater determinants, any slightly modification of the state |$s) can be described in terms 
of a corresponding infinitesimal change of the single-particle density matrix elements. Let us indicate such a change 
by 8p a p. At this point the next task consists in trying to evaluate the change of the total HF energy in terms of these 
quantities. 

Since |$s) is a solution of the HF equations, 14 the first order variation in the energy must vanish so that the 
problem at hand is reduced to determining the sign of the energy change to second order in the 8p a ps. The relevant 
expression is therefore given byJ^ 



A^EffF = - ^ Spa/3 13 — SpfSa + - ^ (v a/3l S ~ V a pSy) 5 p a s5 pgj , 
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where the notation (a, f3) means that only states situated on opposite sides of the Fermi sea are considered in the 
summation. In this formula, we indicate the Hartree-Fock eigenvalues as e a and the corresponding occupation numbers 
as tLq . 

The next step in our procedure consists in constructing a Slater determinant for which the HF energy is lower than 
($>s\H In order to do, we follow Overhauser's idea and choose the new one-particle states to be suitable linear 
combinations of chiral plane waves states situated near opposite points on the Fermi surface, the distortion being 
limited to a very narrow strip. The width of this strip will play the role of a variational parameter. We then carefully 
devise an expression for the wave vector dependent amplitude of the coupling between the plane waves and construct 
the corresponding Slater determinant. In the last step, we calculate the change of the Hartree-Fock energy due to 
this perturbation to leading order in the distortion amplitude from Eq. (1121) . 




FIG. 2: Case of unit chirality when the states of the lower chiral subband are occupied up to the lower threshold of the upper 
subband. The arrows indicate the states that are coupled by the distortion along the k x axis. 



A. Instability for the case of chirality equal one 

Because it presents a formally simpler problem, the first case to be treated is that in which only the lower subband 
is occupied while the upper one is just about to be filled. 15 This situation is depicted in Fig. [2j 

Let us build the new trial wavefunctions as mixtures of the wavefunctions corresponding to wave vector k with 
those corresponding to wave vector k ± Q, i.e. 



^k+QM^k+Q,- 



(13) 



Here, as anticipated, \Q\ = 2kp~ as shown in Fig. [2] 

In evaluating the HF energy change, only the states situated on opposite sides of the Fermi sea are relevant. We 



will consider 



k±Q 



> kp and k < kp. Here, the occupation numbers are = 9 (kp — k), while the amplitude 



satisfies the condition: A^ ± ^ g = The only non-zero matrix elements of dp have the form: 



^Pk4-S. **P%-5. k+G 



(14) 



These wavefunctions indeed describe SDW/CDW-like states. A simple calculation shows that retaining only the 
linear order in the amplitude of the distortion, the spin and the charge densities exhibit spatial oscillations with wave 
vector 2k p. Specifically: 
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SS z (r) 



1 — cos 



0^ 



cos U ■ r 



(15) 



The change in the HF energy is obtained by substituting the expression of the non-zero density matrix elements 
from Eq. (|14l) into Eq. (TT2"1) . The resulting expression can be expressed as: 



A^E HF = A^E HF 



A^E 



H ^HF 



where we have defined 



A^E HF = 

k 

***** = !e6m.+.m 

k,p 

x K+|+^_|)(v# 



+ ,p+%, + ,k-%,+ 



(16) 



(17) 



v k- 1 , +iP - + 4, +iP -_Q, + ,fc + 4, + ) 



(2)£ HF 



A:, /J 



(18) 



(19) 



The Hartree and exchange terms in Eqs. (|18[) - (fH))) contain combinations of the matrix elements of the electron- 
electron interaction. By employing Eq. (J8)), after simple algebraic manipulations, we obtain: 



and 



2L 2 



1 + COS (<^ +| - + COS (<^ + | - ^_|) + COS («^ +| 



(20) 



,p-_f, + ,jfc-f, + Z?+f ,+ + V k-3., + ,p+%, + .k+%, + ,P~%, 



k-p\ 



2L 2 



1 + cos - + cos (<^_g - -) + cos (<^ +# - + - ^_|)] . (21) 



In order to explicitly evaluate the change in the Hartree-Fock energy, we need to assume a specific expression for 
the distortion amplitudes. As a first condition, we will perturb only a narrow region near the Fermi surface. Following 
the same pattern of the proof of reference^, we propose for the present problem the following educated variational 
guess: 



0<4) ! 



n ^ —n 



— T 



^#3, bk F<k <,bk F 

rem c/i- \ r 1 



(22) 



0, 



otherwise 



where b <C 1 is our small parameter and the second arbitrary > 1 . This expression is intentionally chosen to display 
singularities for k — and cf>k — 0. These singularities are crucial to the present proof. 

We now notice that in the expressions of the interaction matrix elements, there appear factors of the type: 
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/ \ (pcosd)ff + k F ) (k cos 6r + kp) +pk sin d>0 sin <hr 
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2k F k cos (j>k 



(23) 
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Since we are only interested in the leading order expansion with respect to b, these cosines can be simply taken to be 
equal to unity, since in the region where the amplitude of the distortion is non- vanishing, both p/kp and k/kp are of 
order b. Accordingly we will assume 




At this point, we recall that k + Q and k — Q must lie on opposite sides of the Fermi sea (i.e \k — ^| < kp and 
+ > kp ), which implies that |cos^>g| > |. 

We can now proceed to the evaluation of the three components of A^Ehf from Eqs. (fT8|) - (fl9|) to leading order in 

b. 

For Aq Ehf the first step is to calculate e? q ~ e t q ■ Using fj 10[) . with energies expressed in Ry, x = i~ and 
x' — we have: 



4zr Aau 1 



l\ rim 1 AAJ 



x' (1 + cosv?') dx 



r s Trr s J J + x' 2 - 2ux' cos u>' 

v 



4u 2 4aw 1 



x' (1 + cos <p') dx'dcj)' 



s nr s J J , / u 2 _|_ x a _ 2 UX > cos 

n n v ' 







u—^J\+x 2 +2x cos <fi^ 



(25) 



The quadratures appearing in this expression can then be manipulated by making use of the results of Appendix [A] 
The result is: 



16a; cos < 



1- — - — In £ccos< 

Z Z7T 



(26) 



where the logarithmic term accounts for the divergence of the derivative of the HF single particle energy near the 
Fermi level. Here, £ is a constant approximately equal to 0.51. 
By substituting (gHJ) and ^1 in dTTJ) we obtain: 



A { 2) E HF = — " / — / d^mr. ( i_ Zli _ -2-lnl^coswl ) , (27) 
(In!) r 2 7r 2 i != 7 cos^ V 2 2tt 15 ^ ' V > 

\ b) s b 

an expression that, to leading order in 6, reduces to: 

where A is the number of particles. 

The Hartree term A H Ehf (containing vq) can in turn be evaluated as follows. By making use of the assumed 
amplitude in (|22|) . we write: 



16b 3 L 2 kp f dx f simp 



6s 6? 



A$E HF = imb \ \*L\% f dlfix ^f I d^^S. (29) 
7r 2 (lnl) rJ VxJ VV J cos<j) x J costp y 

V b > s b 6 

At this point, using the result (|F35|) . the leading order in b of this quantity is given by: 

aS » Ehf k "'y-"' . (30) 
-K 2 r s 
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The last term of (1161) . the exchange energy contribution, is clearly negative and therefore will certainly lower the 
energy. Its evaluation is formidable, for it involves several complicated and seemingly difficult quadratures. Rather 
than attempting to actually calculate it, we will establish a lower limit for its magnitude. 

We will restrict ourselves to the region in which both angles 4>? and <j>p are in the first quadrant. Since this excludes 
some contributions of the same sign, the exchange energy will be underestimated. We therefore have: 



be 6? 



> 



2Nb 3 f ' f ' f f 1 1 Vising 1 y/\sm ip y \ 



dx dy / dip x / dify 



TT 2 (lnl) 2 r s J J J J y/x 2 + y 2 ~ 2xycos (<p x - <p v ) |cos^| y/y |cosp v | 

' b b 

(31) 

It is simple to see that this expression will turn out to be proportional to (In b) . This is due to the presence of three 
singularities in the denominator of the integrand. Of these one stems from the divergence of the Coulomb potential, 
while the other two come from the upper limit of the angular integrations. 

Another simplification is provided by the use of the inequality: 

1 1 

> ,. ■ (32) 



^Jx 2 + y 2 - 2xy cos {<p x - <py) yjx 2 + y 2 - 2xy sin ip x sin ip y 

Using the same changes of variable as in Appendix IB! i.e. t x = t&ii(<ft x /2) and t y = tan(<p y /2), we can rewrite this 
integral as follows: 



> 



l6Nb 3 



6? 



(In I) 



dx 



bq 



dy 



I dt x 



dt„ 



8xyt x t y 



l -tl\ 1 + ^1 



t 2 ' 



l+t y - 



It is clear now that the main contribution to the integral comes from the region around the upper limit of integration 
for both t x and t y . In order to retain the leading order term, a good approximation will be to replace both t's with 
1 — 6/2 in the denominator of the first square root. In this way, we can separate the angular integrations in (|31|) to 
obtain: 



> 



2Nb 3 



bq 



7T 2 (In • 



dx 
fx 



b S 



dy 



1 



W j x 2 



y 2 2xy (1 - f ; 

The last two integrals are evaluated using (|B5I) and other changes of variables (x = ub, y = vb) : 



df, 



dt, 



1 



f .1 



1-ti 



1 + tl 1 - t 2 y V 1 + tl 



(34) 
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8Nb 3 
71 r s 



du / dv 



2 b 2 



(35) 



l (u 2 -v 2 Y + ^- 

The last quadrature is calculated in (|B8I) . leading us to a very simple inequality for the exchange contribution: 



8Nb 3 , , 1 
> ^5 — In ^ In - 

ir z r s b 



(36) 



This term contains a logarithmic factor In i , which allows the negative change in the exchange contribution to control 
all the remaining terms. This concludes the proof for this case. 

The same chain of arguments does apply to the case in which the generalized chirality is greater than one. The 
coupling that produces this kind of instability is schematically shown in Fig. [3J All the formulas we derived in the 
previous case do still apply. The only difference lies in the lower integration limits of Eq. (|25|) . but no relevant 
contribution ensues from this. The matrix elements related to the Hartree and exchange contributions are the same, 
and, as a consequence, the leading order approximation is identical. 
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FIG. 3: HF instability: symmetry breaking coupling for the case in which only the lower subband is occupied. 

B. Instability for the case of chirality less than one 

The argument of the previous Section can be applied when the chirality is less than one, i.e. when both chiral 
subbands are occupied. We can try to break the symmetry by coupling states with the same chirality as well as states 
with different chiralities. When same chirality states are coupled, there is nothing new, as one simply just adds a 
chirality index to the various quantities. In this case, the wave vectors characterizing the oscillations are given by: 
Qft = 2k p and the trial wavefuntions can be written as: 



^kfj, — ^kfi + ^k+Q^S^k+Qf.n + ^k-Q^.kfi^k-Q^fj, ' $7) 

This type of coupling is depicted in Fig. |U 

The corresponding distortion of the components of the spin density and the number density can be again calculated 
up to the first order in the amplitudes: 



SS X (f) 
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k . f.i 
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an Q^-f+J2 A k„ ( sin <^ + ^ 



sin<z>_ 



fc-- 



sinQ, 



E%, 



k.f_t 



cos^ +% 



COS 0, 

k- 



cos ■ r 



cos Qn ■ r 



k . / 1 



fc-- 



4E^(i- 

k . ft 



COS 



cosQp 



cos • r 



,(38) 

,(39) 

(40) 
(41) 



We proceed in this case by choosing an amplitude not unlike the one assumed above: 



A %ti 



o. 



fc3 



\/| sin <f>k\ 



, bk F <k< qbk" t 
otherwise . 



(42) 



The proof of the corresponding instability theorem proceeds then in exactly the same manner. 

As anticipated, there is also not much difference when we try to couple states with opposite chirality (see Fig. [5j. 
Although some of the expressions involved in the derivation do change, the main features of the argument remain 
unchanged. The coupling vector in this case is given by Q = kp + kp. Here, we try to find a lower energy state 

by coupling wavefunctions with wave vector k with those with wave vector k ± Q and opposite chirality. The trial 
wavefunctions then read: 



fe+ — ^fe+ + ^k+Q^k+Q - + ^k-Q^^k-Q,- 



(43) 
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For Aj: ± g = the only non zero variations of the matrix elements of the single-particle density matrix 

operator acquire the following form: 



In this case, the new state is characterized by a similar spin and density modulation: 



5S x {r) = fiY. A k [cos(Q-r 

k 

dS y (r) = h^Aj: [sin (-Q ■ ; 

k 

SS t (f) = hJ^A^ [sin(g-r- 

k 

Sn(r) = 2j2A^\sm(Q-f f - 



J h4-<3 tl Q 



k— ^ 



sin Q ■ r 



sin (Q-f-^ + | 



The corresponding terms in the Hartree-Fock energy change are: 



(44) 

(45) 
(46) 
(47) 
(48) 
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A^E hf 



i i w • "r o / 

fc + — >/* k--2,-nJ n 



6 k-S.fj, C k+S.-fj, 



— n-. 



A 2 

A k j 



I ?V o + ?V o )("-o +«-^o ) Ar: A$ v 



-^,-v,k-^,-n,p+^,v ' "k-%,-n,p+^,u,k+^-pi,p-^,-u 



.p k—-4-,—fj,j 

The relevant interaction matrix elements which appear in the expression of the Hartree term become: 

2wq 



while those determining the exchange contribution reduce to: 



fe+4 



J \k-p\ 

L 2 



(cos (<^ + | - + COS (^_| - + 



fc-H 



COS I 



(^+l-^-#)- cos (^ 

+ cos - + - + i) 

As in the previous cases, we assume the following coupling amplitude: 



At 

k 



In- 



T~ 



\Zl sin ^g| 



o, 



6fc^ < k < ^6fc+, | cos^| < | 
otherwise . 



(49) 



(50) 



(51) 



(52) 



(53) 



(54) 



The first term in the Hartree-Fock energy (fTTl) has a positive contribution which is proportional to b 3 . This originates 
from the same logarithmic factor associated with the divergence of the derivative of the single-particle energy at the 
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FIG. 4: Symmetry breaking coupling of states with the same chirality. 



E 




FIG. 5: States with opposite chiralities coupled to obtain a HF instability 

Fermi level. The Hartree term introduces higher order terms in b due to the presence of the sine factors in its matrix 
elements. Finally, the leading order contribution to the exchange matrix elements is 4ui£ _j. By evaluating integrals 

similar to those in Eq. pip we again obtain a negative energy change of order b 3 ln(l/6). 

V. CONCLUSIONS 

We have been able to formally construct a number of distorted chiral states which, irrespective of the electron 
density, have, within mean held, a lower energy than the spatially homogeneous paramagnetic chiral HF ground 
state, thereby affording a rigorous proof of a generalization of Overhauser's Hartree-Fock instability theorem to a two 
dimensional electron liquid in the presence of linear Rashba-Dresselhaus spin-orbit coupling. It is important to notice 
that, as mentioned in Section [ill! to establish the instability we have not allowed for momentum space repopulation: 
inclusion of this phenomenon would have further lowered the energy of the trial states while greatly increasing the 
difficulty of the analysis. 

It is worth mentioning that the states that have been analyzed in this paper differ in a number of ways form the 
original spin/charge density waves proposed by Overhauser. Our states are chiral density states, that display both 
spin and charge modulations. The presence of charge modulations cannot be ignored (as commonly done in the case 
without spin-orbit) and the Hartree term has to be evaluated explicitly. The exchange gain is shown to be larger than 
the kinetic energy plus the ensuing Hartree terms. In this respect the difference with the plain electron liquid is quite 
marked for it is precisely the electrostatic term that is believed to be fatal to the charge density waves in that case. 
We prove that this is not the case for the chiral waves. The main reason is that the exchange energy behaves quite 
differently as compared to that in the absence of spin-orbit interaction. In particular it does not necessarily favor a 
homogeneous chiral instability as shown in RefAi. Our calculations show that in a vast class of inhomogeneous states, 
the exchange energy can induce a chiral instability for all densities of the electron liquid. Of particular importance is 
that the chiral states coupled by the type of mean helds explored here do always have opposite spins. The fact that 
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such couplings are possible for all values of the spin-orbit coupling is another consequence of the non analyticity of 
the properties of the electron liquid on this parameter. 

Our results only represent the first step in understanding non homogeneous states in this interesting many-body 
system. How the inclusion of correlations will modify the HF scenario is of course a most important question. Advances 
in this respect can in principle be pursued by following the program outlined in Reference [iq. Such a study will require 
a vastly more complicated analysis. An alternative route is to make use of perturbative techniques to establish the 
role of correlation effects in the high density limit. Part of this program has been carried out in Reference [l?] for 
homogeneous states. 
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Appendix A: Elliptic integral expansion 



For the purpose of our calculations the following expression must be evaluated in the limit of small u: 



2ir 1 



x' (1 + cos (p) dx'df 




2tt 1 



•y/z 2 + x' 2 — 2zx' cos if 
x' (1 + cos if) dx'dif 



z — l+u 



\J z 2 + x' 2 — 2zx' cos if 



An obvious problem is the presence of singularities in the integrand for z = 1. 
Let us define: 



(Al) 



*(-)-/ / fQ- + «*<P)&*P (A2) 

J J \J z 2 + x' 2 — 2zx' cos if 

v 

Of course, A is not differentiable in z = 1. Still we have to evaluate ® A q^ an d expand it in an asymptotic series 
around z = 1. 



dz 



A(z) 



2tt 1/z 

f f y (1 + cos <f ) dydf 
dz \ J J ^/1 + y 2 - 2y COS f> 



2tt 1/z 



y (1 + cos if) dydf 
\Jl + y 2 — 2ycos</3 



/ 2ir 1/z 



dz 



\o 



y (1 + cos if) dydf 
a/1 + y 2 — 2ycosf 



(A3) 



The derivative in the second term can be evaluated in the following way: 



d_ 

dz 



2ir 1/z 







y (1 + cos f) dydf 
\J\ + y 2 — 2ycos f 
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2tt 

1 f y (1 + cost/?) dip 



z 2 







\J\ + y 2 — 2ycos< / 9 



2tt 



.2 



1 /" ( 1 + cos f) dip 

~~2 



(I 



%Jl + z 2 — 2zco%<p 

We can integrate this expression over its angular part so that the result is expressed as: 

(l + cos<^)d<^ 1\z-\\ ( -4z \ 

E 



yjl + Z 2 - 2ZC0S( < 5 2 \(Z-1) 



2 



+ 2 (£+1)!^ f 



^ l*-i| V^- 1 )' 



where J5 and are the elliptic integrals of first and second kind defined as in 18,19 . 
The asymptotic expansion of the elliptic integrals around z = 1 leads to: 



where p — J J -^4= 

n n vi 



2tt 1 

y(l+cos tp)dydtp 



, , _ 5.6639. 

o o v 1 +y 2 - 2 y cas v 

Integrating over z we finally obtain: 



with £ ~ 0.51. 



Appendix B: Evaluation of useful quadratures 

We begin by calculating here the leading order term in the expansion of the expression 



arccos : 







Expanding the upper limit and setting t = tan ^, the integral becomes: 



I ~ 2^2 



1 2 

dt 



i-t 2 y 1 + t 2 





- 1 - 







(A4) 



(A5) 



— ^ = p- 121n2 + 4 + 41n|z- 1| , (A6) 



4(1 + u) - A(l - u) ~ 2 (p- 12 In 2 + 41n|u|)u = 

~8u(ln|£u|) , (A7) 



/= / ^El^ forb«l, (Bl) 

' COS</7 



^ ' — + TTT>' (B2) 



The divergence in the limit 6 — > stems only from the first term and we therefore proceed to try isolating the 
singularity: 
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(B3) 



V2 / dtJj+p (In (l-t))' + 0.5256 + 0(b) . 
o v 



An integration by parts is used in the remaining integral to further isolate the singular contribution: 



21 On(l-t)) 



l + t 2 



+ V2 y dt^yZ^Z^ in (l-t) + 0.5256 + 0(6) . 



(B4) 



At this point the non singular second term is calculated numerically so that the final result is: 



J~ In 6 + 0.9475 . (B5) 
The following integral is used in the proof of the HF instability theorem: 

J=ff , ^ (B6) 



for b -C 1. 

Because the singular behavior originates from the region where u ~ u, in order to find the leading order term, we 
approximate u 2 — v 2 with 2u(u — v). Since the function is symmetric with respect to the interchange of u and v we 

can use the relation J du J dv f(u,v) — 2 J du J dv f(u,v). 



Then: 



J ~ 



v 7 ? u vf 

du f dv 



{ U { J(u-vf + ^~ { u V2V2 



/du ^ ( bu 
~ n V2Vl 



V? 



u 



+ l-ln u-l + J(u-lf + ^). (B7) 



l 



The integrand of the first term has a logarithmic singularity in the limit of small b while the second one is instead 
analytic in this limit. The main contribution to the integral is therefore: 



' bu \ 
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